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Use of a cDNA microarray to 
analyse gene expression 
patterns in human cancer 

Joseph DeRisi'*, Lolita Pcnland- & 
Patrick O. Brown- (Group I); 
Michael L Bittncr^*, Paul S. Melizer\ 
Michael Rax-^. Yidong Chen\ Van A. Su^ & 
Jeffrey M. Trent-^ (Group 2) 



The development and progression of cancer^'^ and 
the experimental reversal of tumorigenicit/-* are 
accompanied by complex changes in patterns of 
gene expression. Microarrays of cDNA provide a 
powerful tool for studying these complex phenom- 
ena^®. The tumorigenic properties of a human 
melanoma cell line, UACC-903, can be suppressed 
by introduction of a normal human chromosome 6, 
resulting in a reduction of grovy^h rate, restoration 
of contact inhibition, and suppression of both soft 
agar clonogenicity and tumorigenicity in nude 
mice^-^*®. We used a high density microarray of 
1,161 DNA elements to search for differences in 
gene expression associated with tumour suppres* 
sion in this system. Fluorescent probes for 
hybridization were derived from two sources of eel* 
lular mRNA (UACC-903 and UACC-903(+6)] which 
were labelled with different fluors to provide a direct 
and internally controlled comparison of the mRNA 
levels corresponding to each an-ayed gene. The flu- 
orescence signals representing hybridization to 
each an-ayed gene were analysed to determine the 
relative abundance in the two samples of mRNAs 
corresponding to each gene. Previously unrecog- 
nized alterations in the expression of specific genes 
provide leads for further investigation of the genet- 
ic basis of the tumorigenic phenotype of these cells. 

DNA microarr;iys« contain ini; IJftl total clcmenis. 
including S70 dilTcrcnl cDNAsaiul controls'^*" (sec 
Methods), were printed roboiically onto a jthiss micro- 
scope slide in four quadrants coverini; an ;irea of about 
1 cm- (Fig. I ). We prepared fluorescent cDNA probes 
using total poly (A)* mRNA from UACC-9()3 cells and 
UACC-903(+6) cells by labelling with a green and red 
fluor» respectively. A mixture of the lw«> flourescently 
labelled probes wa.s hvbridi/.ed xo the DNA microarrav. 
This comparative hybridization method, coupled with 
the doping of synthetic standards and an estimation of 
statistically significant deviation for local background 
variance allowed a direct and quantitative comparison 
of" the relative abundance of individual DNA sequences 
in this complex Siunple'*"**. We added a set of synthetic 
poly ( A)*-lailed niRNAs' to the purified mRNA from 
each cell line as internal standards to jtssisi in quantita- 
tion and estimation of experimental variation intro- 
duced during labelling and reading. Targets 
complementary to these standards were included, in 
duplicate, on the micronrray. Based <mi these standards. 
mRNA species comprising 1 : lO.OOO of the mass of the 
poly (A)* RNA could readily be detected. 

In a representative two-coU)ur fluorescent scan of all 

I ir.i . . I «...,.... ^r: . .\ .1. . ... . t, ... 



spend to genes preferentially expressed in the tumori- 
genic UACC-903 cell line, and the reddish spots corre- 
spond to genes preferentially expressed in the 
non-tumorigenic UACC-903(+6) cell line. Genes 
expressed at approximately equal levels in the two cell 
lines appear yellow or brown. A portion of the array at 
higher magnification highlights the diverse pattern of 
differential expression observed (Fig. 26). In Fig. 2c, rec- 
tangles corresponding to specific array elements are 
coloured to reproduce the hue and intensity of the fluo- 
rescent signal at each element. The hybridization signals 
from a duplicated set of genes are shown juxtaposed, to 
illustrate the reproducibility of the hybridization signals 
for each gene. 

To address the possibility that an apparent difference 
in expression might result from experimental variables 
unrelated to the difference in chromosomal composi- 
tion between the two cell lines, we examined the vari- 
ance in expression for 90 ^housekeeping' genes. We 
selected the.se genes ba.sed on the assumption that they 
would not be differentially expressed between the two 
cell lines. The averaged red/green ratio for this subset of 
genes was 1. 1 3. The averaged red/green ratio for the set 
of five internal standards was 0.97 ( 11 = 10). The vari- 
ability in the expression level of the housekeeping genes 
probably overestimates the experimental variability in 
measuring differential expression. As a conservative stan- 
dard, an absolute fluorescent signal (red or green) with 
an intensity greater than that observed at the contr 1 
array elements containing total human genomic DNA 
was considered to represent specific hybridist ion. Gene- 
speciilc hybridization was therefore only considered sig- 
nificantly different between s;unples if the following tw 
criteria were met: i) the signal intensity (green or red) 
exceeded this threshold; and ii) the logarithm of the 
red/green fluore.scence signal ratio differed by >3 S.D. 
from the mean logarithm of this raiit> for the * house- 
keeping gene panel (that is. ratios <0..S2 or >2.4). 

By ihe.se criteria, mRNA levels for 1 5/870 ( 1 .7%) genes 
were significantly diminished, while the mRNA levels 
for 63/870 (7.3'>o) genes were significantly increased in 
association with .suppre.ssion of tumorigenicity by intro- 
duction of chromosome 6. To test the reliability f 
microarray hybridization results in identifying differen- 
tially expressed genes, we analy.sed 16 genes by north- 
ern analysis. In each case, the results of northern analysis 
corroborated the differentia! gene expression identified 
by microarray hybridization (Fig. 3). 

Significant differences in expression between these 
two cell lines identified several genes as candidates f r 
determining features of the tumorigenic phenotype f 
the melanoma cells. For example, among the genes 
detected with significantly higher expression (>lO-fold) 
in the tumorigenic cells was the human brown locus pro- 
tein (TRPI/melanonia antigen gp75). This is the most 
abundant glycoprotein in mclanocytic cells and a critical 
melano.some membrane protein'- ' \ Additionally* its 
expression is reduced when melanoma cell lines are 
induced to differentiate by treatment with HMBA*^''^ 
Also expressed at a significantly higher level was a spliced 
variant nf the mRNA encoding myelin IM.IV1)M20. This 
is widely expressed in neural crest derived cells in early 
development and has been suggested to play a role in 
cell-ceil signaling during develt»pment''*. 
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Rg. 1 Propenies of cDNA microarrays. a. A fluorescem scan of DMA printed onto a poty-lysine coated slide. The DMA is stained with a DNA-speciftc fluorescent 
dye, YOYO. The center-to-center spacing of adjacent spots is 450 m- allowing the potential for up to 10.000 spots/2.54 X 7.62 cm microscope slide, b. Effi- 
cient blocking of hybridization to DNA repeals. Hybridization of ftuorescein-tabetled poly (dT)* to arrays in the absence of competitor produces strong 
hybridization to immobilized poly (dA)* as well as to some cDNAs. such as the EST T64827 shown. Rhodamine-labelted cDNA (red) from the UACC-903 cell 
line hybhdtzed in the presence of poly (dA)' blocker shows little if any signal at either site (Total H = total human). Smitarty. hybndization with fiuorescein-labetted 
Coti DNA in the absence of competitor produces bright signal on immobilized Coti DNA. total human DNA and at some cDNA elements (presumed to con- 
tain highly repeated sequences, such as R23416): while Rhodamine-tabeiled cONA (red) from the UACC-903 cell line produces little it any signal at these 
locations when hybridized tn the presence of excess unlabelled poly (dA)*. and human Coti DNA. The absence of signal at some cDNA locations following 
UACC-903 cDNA hybridizations also indicates that the PCR-amplified. plasmid vector sequences at all cDNA targets do not contribute signiticant hybridiza- 
tion signal, c. Schematic of the array organisation. Robotic printing from 96 well microliter trays was earned out with 4 print heads, spaced to fit into 4 adja- 
cent microliter wells. This maps the contents of each tray into four separate quadrants on the glass slide. A colour-coded map of the general distribution of 
target types in each of the resulting quadrants is shown. 



cLs were clcv;iiai by the ;iddiluni of a norniiil chromo- 
mmc 6 ( 17 jicncs) arc kmnvn lo be activated by IFN-y, a 
cardinal proinflammatory cytokine that, among other 
activities, induces expres.sion of I he gene products of the 
M HC classM locus. For example, theniRNA enci>ding_ 
monocyte chemotaciic protein I (MCAF/MCFl ), a 
cytokine that induces mont>cyte chemotaxis and activa- 
tion'^"', was more than 10- fold le.ss abundant in the 
lumorigenic cell line. In the .skin, MPCl is critical in the 
regulation of cutaneous monocyte trafficking'*'"'', and 
elevated expression plays a role in suppression of tumour 
growth and metastasis'***-'. The mechanism by which 
these interferon-y regulated genes are induced in UACC- 
903 cells by transfer of a normal chromosome 6 remains 
to be determined. It is worth noting, however, that the 
interferon-7 receptor gene is localized to the distal long 
arm of human chromosome 6. 

Finally, several genes that showed > 10- fold higher 
expression in the suppressed UACC-903 (-♦■6) cells have 
previously been recognized in other models of tumour 
suppression. Ni st notably, there was elevated expres- 
sion of the mRNA encoding WAFI (p2l ), a key media- 
t r f tumour suppression by p53 (ref. 18). The p21 
protein had previously been identified as a melanoma 
differentiation-associated antigen (termed mda-6)'^*-*'. 
In melan ma cell lines suppressed for metastasis by the 
introduction of chromosome 6, expressi n of WAFI 
(p2l) mRNA and protein correlates inversely with 
metastatic potential-". 



These results provide a wide view of the diverse sys- 
tems that arc altered in this model system of tumori- 
genicity» and focus attention on specific gene products 
and pathways that may be of particular importance in 
this timiour type, _ _____ 

Our ability to classify human cancers in a way that 
reflects the underlying molecular pathology or that 
anticipates their potential for progression or response 
to treatment, remains primitive. Using cDNA microar- 
rays to dcfme alterations in gene expression associated 
with a specific cancer may be an efficient way to uncov- 
er clues to the specific molecular derangements that con- 
tribute to its pathogenesis and thus identif)' potential 
targets for therapeutic intervention. Moreover, recogni- 
tion of pathognomonic alterations in gene expression 
might provide a basis for improved diagnosis and mol- 
ecular classification of cancers and thus allow selection of 
the most appropriate therapeutic strategies. 

Public databases of human expressed gene sequences 
contain partial sequences of at least 40,000 different 
human genes", and efforts to develop a human tran- 
script map have developed rapidly-'. Based on the high 
yield of information obtained using an array of < 1,000 
different genes, a more comprehensive survey of gene 
expression patterns, using a more complete array of 
human genes, will likely provide a rich source of new 
and useful insights into human biology and a deeper 
understanding of the gene pathways involved in the 
pathogenesis of cancer and other diseases. 
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Fig. 2 DNA microarray analysis of changes in gene expression between the tumorigenic cell line, UACC-903. and its non-tumorigenic denvative UACC903(*© 
f^"''f?.tV^^L^^^^^T °^ ^ c^iromosome 6. a. A ratio image of the results of simultaneous hybridization of Rhodaminel lO-tabelled cDIMA (qreen) 

from UACC.903 and CyS-tabelied cDNA (orange-red) from UACC-903(*6) to a microarray. To produce this .mage, the scan images corresponding to eachX 
orescent probe were combined as the appropriate colour channels m a single image. Arrows indicate the location wttn.n the array of the corresoondino 
genes analysed by northern blotting (F.g. 3). b. A magnified image of the area of the array boxed m white m (a), c. Maomfied image of three cDNAs indentiV 
^ed by arrows inja). representing the cDNAs for: left. MCAF/MCP- 1 (r/g ratio >10): centre, fi^act.n (r/g ratio 1 .04): and^.ght, »- 7 .anf/chymof/yps/n (r/g ratio 
0.2) (see Fig. 3]. d. simplified representation of ratio hybridization results. Quanhtaiive fluorescence intensity data is extracted from each anay taroet The aver- 
age target colour ratio determines the hue of each box and the average intensity determines the brightness of each box. tn this image the ordeTof'the boxes 
corresponds to their ongmal order in the microliter plate from which they were printed. Duplicate printings of the same plate can be 'examined side bv side 
as in the first two rows shown here, to assess reproducibility of the hybndizaiion results (see text). Numbered arrows indicate the location within the arrav cor^ 
responding to genes analysed by northern blotting in Fig. 3. ' 



Methods 

Generation of microarrays, hybridization, scanning. Ilu* 
prep.u.iliiin (il tiUlcil inicnvsCiipc .sliiif.s ami Mil»Ci)tit'm lohui- 
w printint; ol DNA waN c.irrial »)Ut m a inaniur Ntinilat id tliai 
ilcMTilvit . Hricllv. pre-(.kMiK*il gl.ts> sliiks wm ircaicil uiih 
poly. I -lysine Miltilidti tSicnu) In lorni an atllicsivc Mirlaa* lur 
printinj:. I't'li pnulucis, purifial hy ctlianol puril'icalion. uviv 
rcsuNpcmlai in 3\ SSC. A ciiNi»»ni Iniili arrayini; mlmi picktil 
up aiut iicposiicd .small vtiUiinos nannliier.s) nl t)NA «»nio 
the .sillies. Alter printing, the slides uvrc washai in a i),y'/u SI ).S 
.st)UilH»n. The renuiining iuiund DNA was ilenaimeil hy .suh- 
merging the Nliiles in 93 "C Jisiilleii water tor 2 inin Inllnweil hv 
a hriel iva.sh with ethanol. DNA was L'\' eii»sslinkei1 to (he 
slides (Straiageiie Stratalinker. W) ml). To prevent non-specilu 
probe binding, the .slides were blocked by rinsing m a sol m ion 
t>r7(>-niM succinic anhydride dissolved in 0.1 M b<»ric acid pH 
S.n. ciHiiaining 3.=^% I -methyl-2.pvrr«»lidinone (Atdrich). 
Addiiiniial pniU»cols and parts lis! pertaining it) iuicroarra\ 
fabrication can be obtained from hi rp:y/cnigiu. sianford.edu/ 
pbrown. 

Purified, labelled cDNA was resuspended in 1 1 pi i>r3..Sx SSC! 
containing 4 fig of pi)ly ulA)* DNA. :.5 L tv/; tUNA, 4 fig of 
human Cot I DNA ((iibco HKI.l.and 0.3 fd of lOV.. SD.S. l'rK»r to 
hybridi/aiion. ihe soluiion wa> boiled )«>i 2 inin then allowed lo 
coot lo loont leinperaiiire. I lybriili/aii4Mi was carrieil out ai 
h2 \' lor *U h in a water bath. Prior lo scanning, slides were 
washed in > SS( !. 0.2".. SDS Itir 5 min and 0.2v SSC for I inin. 

local l.iMT inurosv.i»pe Iniili bv S. Smith with sofiwaiv wiiiieii 



in N. /i\. A separ.ue scan, using the appro[niate excitation line, 
was done tor eacii ol the two tluoroplH>res used. Data was col- 
kcied .11 .1 in.ixiimini ivm»IuIiom t>l 9 inicioiis/pixel with 12 bits 
ol depth 

Probe preparation and labelling. UNA was extracted from cells 
using the IVia/ol reagent (ITl Inc.). fc»Ili»wing the manufactur- 
er's direct ions. cDNA piiibes were synthesi/ed fn»in singly oligo 
dl-selecied {Phainiacia) mK\A pools. Huorescenlly labelled 
cDNA was prepared from niUNA in «»ligo dT-priined polynier- 
i/aiion using SuperScnpi II reverse transcriptase (LTI Inc.). 
I he poi»l of nuclei) I ides in the labelling re.iciit»n was 0.5 mM 
dC.TP. dAlP and dc:TP and 0.2 niM dTTP. Phiorescent 
nucieoiides. Kh.ulamine I lo dUTP (Perkin F.lmer Ceius) or 
Cv} dL*TP ( AnuTshamI, were present at 0.1 inM. Probes were 
purincd bv gel chromaiogr.iphy (llioSi»in ft/KioRad) and 
eihanni piecipttaiion. 

Selection of cDNA elements and generation of control Icm- 
plates. Synihetic cDNAs were prepared by cloning random 
/iiii//HI and Hi//dlll ended fragments of coli DNA in the vec- 
tor pSPM poly (A)* (Proinega). linearizing isolated plasmid 
DNA \Mi\\ KtvlU and synthest/mg poly (A)* tailed UNA COm- 
pleinentary to the insert hmn the resilient SP(i promoter 
(PoMiiega). Prior In use. the s\ntlK'si/ed KNAs were .selected on 
tihgo d I' cellulose. The largext group of cl )NAs consi.steil of 674 
cD.V-V clones |i,mi the iMli .irraveil norin,iIi/ed infant brain 

l..;j.ii ; . I hc-.>c (.t.-;;. . . .... sv.^^i^J i,. i; ;er\- NIB 

hbrarv member ih.it conespoiuleil tn a named gene according 
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Fig. 3 Nonnern hyoridization suostantiat* 
ing the consistency ot the cDNA microar- 
ray resutls. Corresponding locations within 
the cDNA microarray illustrated m Fig. 2a 
are provided lor 1} Waf-yp21: 2) MARCKS: 
3) co//ageaase: 4) MCAF/MCP-l; 5) f^T- 
anfichymofryps/n; and 6) ^ -actin. The sig- 
nal detected by a radio-labelled ^acttn 
probe represents a control for loading vari- 
ance, with a red/green ratio observed on 
the cDIMA microarray (Fig. 2a,c) for fi-actin 
ot1.04. 

to the UniGcnc EST cliiMcring sys- 
tem-*-*. The sccnnii largest grmip o( 
iTloiU'S ctin.sisicil (»t* 1 83 wqucnccd cI^NA 
clones i;cneraiec1 In Mibiraciinn of cDNA 
Irom ihe chronn»Mimc*fr supprcsNed 
non-iumoriitenic UACCW(I3 (+6) cell 
line with cDNA frc»m its parental liimori- 
genie cell line UACtC-yo.A (ref. 9). 
Approximaiely I (Ml aiiitttional genes 
(u»tal X7() genes arrayeil) were ohtaineif 
iVoni HST libraries on the kisis of their 
a- 1 'Antichymotrypsin expression paiiern ( lissiie specific, anil so 
ttn). luich array inchiiietl the lollowing 
hyhriJi/aiion conlrols: plasmiO vecitir, 
laniKta. OX I7-) phage, total human DNA. 
human C'.oll |)NA. ami |>c»ly {AW The 
syniheiic siamlanls usi\l lor normali/.i- 
lion ol signals in each wavelength were 
also arrayet). ('ontrols were incluileil in 
each ijiLuirani of the array to assess ihe reproilucibiliiy oi the 
hvhritii/aiion signal. Twt» plates til'cDNA clt»nes (ileriveil iVoin 
the UAC;(;-*)il3 >uhiracteil library) were alsii arrayeJ in ilupli- 
cate. riileliiy tit' the I'nigene array relative i»» tIbKST was lesieil 
by set]uencing ot a raiulom s;imple i»l II cli>nes usetl U>r 
microarray const ructiitn. All sei^ueiices were iileniical with ihe 



MARCKS 



Collagenase 



MCAF/MCP-1 



ip-Actin 



correst'Niiuling Jbl-..ST entries. AilJilionally. ea<.h mum.ni.i\a! 
cDNA Irom the l^\Ct:-V«i3 subtracieil library %\a^ Nei|neiKe.!. \ 
listing ol cl>NAs comprising this microarray whkh ueic 
ilerivcil I rom I he Unigene and housekeeping' panel ».aii l>e 
obtained from hiip://w'w\v.nih.gov/l>lU/lCC/AKI\.A^ 'C\pn - 
html. 

Northern blot analysis. Total RNA, 10 pg (vr lane, was elec- 
trophnresed in I. J'Siagnrose- formaldehyde gels and transferred 
onto nylon membrane (Hybond-N*. Amcr.sham) by capillary 
blotting overnight. For DNA prolvs insert fragments from the 
Soarcs IMH cDNA library'" were tibtamed by veci»»r PGR for 
p2l, MAUCKS. a-l-amichynu»irypsin and (J-actin. Pu»be> lor 
fibroblast collagenase and MC.AF/NU.P- 1 were isolated from a 
UACC-y03(+(>) enriched cPNA library*' with all probes 
labelled by random priming. Tillers were washed to a >iiin- 
genq* of (I.I x SSC at A2 lor :() min. 

Web sites. hilp://cmgm.stanl\>rd.edu/pbri»wn for pnuocois and 
parts list pertaining to microarray fabrication, 
ht I p:// www. nchgr.nih.gov/ 1 )IK/l.('.(i/AKUAY/expn.html lor a 
listing of cDNAs comprising ihis microarray which weiv 
derived from the Unigene and 'htuisekeeping panel. 
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